IN-PLANE CONDUCTIVITY ANISOTROPY IN UNDERDOPED CUPRATES IN 

THE SPIN-CHARGE GAUGE APPROACH 
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Applying the recently developed spin-charge gauge theory for the pseudogap phase in cuprates, 
we propose a self-consistent explanation of several peculiar features of the far-infrared in-plane AC 
conductivity, including a broad peak as a function of frequency and significant anisotropy at low 
temperatures, along with a similar temperature-dependent in-plane anisotropy of DC conductivity 
in lightly doped cuprates. The anisotropy of the metal-insulator crossover scale is considered to be 
responsible for these phenomena. The obtained results are in good agreement with experiments. 
An explicit proposal is made to further check the theory. 
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The normal state properties of cuprate superconduc- 
tors have shown a number of unexpected features, partic- 
ularly in underdoped regime, where the pseudogap effects 
are pronouncedi One of the most recent surprises is the 
strong in-plane temperature-dependent anisotropy of DC 
conductivity, observed in lightly doped untwinned LSCO 
and YBCO crystals^ The maximal anisotropy pb/pa — 1 
reaches 50% for 3% doped LSCO which is far too big 
compared with the orthorhombicity (up to 1.7%)^ A 
"natural" explanation of this unusual behavior would 
be the "self-organized" charge stripe structure, proposed 
by a number of authorS)2i4which was also suggested to 
be responsible for the occurrence of superconductivity^ 
The conductivity is indeed higher in the stripe direction 
(along the a-axis)£ Meanwhile, Dumm et al^ recently 
measured the in-plane AC conductivity of untwinned 3% 
and 4% doped LSCO crystals. In the far infrared region 
at high temperatures, the data are consistent with a sim- 
ple Drudc model. Below 80 K, a broad peak appears at 
finite frequencies (f2 ~ 100 cm" 1 ) bearing a close resem- 
blance to the peak found in temperature-dependent DC 
conductivity for the same composition; a significant a- 
b resistivity anisotropy is observed in complete analogy 
with the DC c&sem- These AC data seem also to support 
the presence of charge stripes. 

However, there are several substantial difficulties in the 
intuitive "rivers of charge" interpretation: The mean- 
field theory predicts the statically charge-ordered state 
is an insulator j£ while experimentally these lightly doped 
cuprates show metallic behavior at high temperatures. 
Also, the observed anisotropy ratio is too small compared 
with quasi-one-dimensional (Q1D) conductors, usually 
showing order-of-magnitude bigger conductivity in the 
chain direction. To avoid these difficulties, the "elec- 
tronic liquid crystal" scenario of meandering stripes^ 



is invoked to induce metallic conductivity and to reduce 
the expected anisotropy. This does not solve the prob- 
lem, either. In fact, a closer examination of the data2*£ 
reveals that the anisotropy effect is most pronounced in 
the limit lo,T — > 0, in contradiction with the fluctuat- 
ing stripe picture: One would anticipate much bigger ef- 
fect of anisotropy at some characteristic frequencies of 
the stripe fluctuations, rather than the static limit. This 
was also pointed out in Ref. 0. Also, as evident from 
the experimental curves the major source of the in- 
plane conductivity anisotropy is due to the shift of the 
Metal-Insulator Crossover (MIC) scale, and a stronger 
anisotropy is found in the localized, instead of metallic 
regime. Up to now, many authors attribute this MIC to 
the disorder effect^ or more specifically to localization 
in Q1D systems^ However, there is a fundamental diffi- 
culty in that approach: There is only one mobility edge in 
disorder-induced localization for anisotropic systems^, 
at least within the scaling theory. This means the system 
cannot be localized in one direction, while delocalized in 
another. 

Recently, we have developed a spin-charge gauge ap- 
proach to describe the pseudogap phase in cuprate 
superconductors, particularly focusing on the MIC 
phenomena>ii*i£ In this approach^ based on spin-charge 
decomposition applied to the 2D t-J model, the spinon 
dynamics is described by a non-linear cr-modcl with 
a theoretically derived mass gap m s ~ J(S\\nS\) 1 ^ 2 , 
where J is the exchange integral, 6 the doping concen- 
tration; the holon is fermionic with "small" Fermi sur- 
faces (cf ~ tS) (with t as the hopping integral) centered 
around (±7r/2, ±7t/2) in the Brillouin zone and a "Fermi- 
arc" behavior for the spectral weight. Both holons and 
spinons are strongly scattered by gauge fluctuations. As 
an effect of gauge interaction, the spinon mass picks up 
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a dissipativc term: m s — * Mr = (m 2 — icT/x) 1 ^ 2 , where 
X ~ is the diamagnetic susceptibility and c a numer- 
ical constant. This shift in turn introduces a dissipation 
in the spinon-gauge sector, whose behavior dominates the 
low energy physics of the system. The competition be- 
tween the mass gap and the dissipation is responsible for 
the MIC, giving rise to a broad peak in the DC conduc- 
tivity. At low temperatures the antifcrromagnctic (AF) 
correlation length £ ~ mj 1 is the determining scale of the 
problem, leading to localizing behavior, while at higher 
temperatures, the de Broglie wave length At ~ (x/T) 1 ^ 2 
becomes comparable, or even shorter than £, giving rise 
to metallic conductivity. Hence in this approach the MIC 
is mainly due to correlation and AF order, rather than 
disorder effect. 

In this Report we generalize this approach to 
frequency-dependent phenomena and show the AC- 
conductivity exhibits a maximum as function of fre- 
quency in an exact analogy with the DC conductivity 
maximum due to MIC. This outcome is fully understand- 
able: In the presence of an external electromagnetic field 
of frequency f2, that frequency will replace the tempera- 
ture T, playing the role of cutoff parameter, for fl > T 
at low temperatures. Moreover, if we assume the AF 
correlation length £ is anisotropic, we argue that as a 
consequence the same is true for the MIC scale, then the 
in-plane conductivity anisotropy as well as the close par- 
allel between the above two sets of experimental data can 
be interpreted in a self-consistent way. 

Let us now analyze the theory in more details. In a 
gauge approach the physical resistivity is calculated using 
the Ioffe-Larkin addition rule>^ p = p s +Ph, where p s and 
Ph are the spinon and holon resistivities, respectively. In 
order to evaluate the relevant current-current correlation 
functions, one first integrates over holons and spinons, 
and finds that the gauge propagator in the scaling limit, 
uj,q,uj/q — ► 0, has a Reizer singularity 1 ^ for the trans- 



verse component (A T A T )(ui, q) 
where k ~ S is the Landau damping. 

The effect of Reizer singularity on gapless fermions is 
subdominant; it has been analyzed in Ref. and at 
finite T it gives a scattering rate for the holon of order 
T 4 / 3 instead of the usual Fermi liquid result ~ T 2 . To 
include, non-perturbatively, the effect of gauge fluctua- 
tions in the spinon current-current correlation functions, 
we expand the spinon propagator in Feynman paths, as 
justified by the mass gap and we integrate over velocity 
fields in the eikonal approximation^ Being gauge invari- 
ant, the correlator of the spinon current depends only on 
the gauge field strength. In the scaling limit only the 
magnetic components Fij are relevant (see Ref. Il2|) ; the 
corresponding propagator at finite T is given by 



(Fij(x)F r s{0)) = [S ir 5 JS - Si s Sj r ] 



dui 
2^ 



dk |fc|V fe - 
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where x = (x,x°) with x° as the time variable. In the 
presence of an external electric field (the probe in linear 
response theory) with frequency f2, the integration over 
uj should be cutoff at \u>\ ^ A =max(T, d). (We further 
assume, for technical reasons, that Ax° <C 1, as justified 
a posteriori, see Ref. Il2|) . 

Let us first focus on the static response and set A = T. 
As u) < T, we approximate coth — ^f, then the uo- 
intcgration in Q becomes 



dio 



2T 



(2) 



eventually leading to an evaluation of the large-scale 
"magnetic" propagator: — i(T / ^ttx)Qt exp (— Q\ |x| 2 /4), 
where Qt = (Tn/x) 1 ^ 3 is the inverse of an anomalous 
skin depth, identified as the length-scale of gauge fluctu- 
ations. The main effect of the gauge interaction on the 
spinon propagator is a renormalization of the mass term 
in the exponent in the limit xq 3> |x|: 



m s x 



X 2 



T 



Yx Q ^ 2 



x\Qt s xl 

1 ™2 ' 



. ( 3) 

where / and g are regular functions, whose explicit in- 
tegral representations are given in Ref. 0. In the eval- 
uation of the spatial Fourier transform of the current- 
current correlator, the shift eventually leads to a com- 
plex saddle point for \x\Qt with absolute value C ~ 0(1) 
and phase factor e" r ' 4 , which in turn introduces a dissi- 
pation term in the spinon gap: 



M T = (m 2 s -icT/ X ) 1/2 , 



(4) 



where ic — /(Ce J7r / 4 ). The competition between the gap 
term m 2 and the dissipation T/x leads to a MIC upon 
the decrease of temperature, yielding a broad peak in the 
DC conductivity for T ~ X m 2 ~ (t/S)\8lnS\ ~ *j In J|, 
thus shifting to lower temperature upon doping increase. 
More precisely the behavior derived for the DC conduc- 
tivity is given hyi— 



a(T) 



/"(C e W4)|M T | 



1/2 



sin(i argA/ T ), (5) 



where /" means the second derivative. 
We turn now to the AC conductivity at T 



and set 



A = 0. In this case coth ^ is replaced by sgnw and the 



-integration in becomes 



duj 



2n 



M 2 + (^) 



(0) 



Up 
the 



to logarithmic accuracy, one finds for 
magnetic propagator at large scales: 



-i(T/4^)Q2Acxp(-Q2 |^ /4) whcrc Qn = (n K / x )i/s 
and 0<A<l/2, as follows from comparing (JSJ and J3J. 
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Repeating the steps of the DC calculations with this 
parameter A included, we find as the analog of 

m a -» M n = (m 2 s - ic\n/ x ) 1/2 . (7) 
For fl <C 2ms , one easily obtains for the AC conductivity: 

gi( "^( Ar(c e w 4) |M n | ) 1/2sin (r rgM ^- < 8 > 

We see that the behavior of the AC conductivity at T = 
is rather similar to that of the DC conductivity with 
a broad peak corresponding to a MIC, hardening and 
shifting to lower temperature upon doping increase, (see 
Fig. 2}. Although the replacement of T by fl as cut- 
off yields a priori only an order estimation, the presence 
of the factor A (mainly coming from the factor coth ^) 
suggests that the position of the peak in the AC conduc- 
tivity is shifted by a factor A -1 to higher frequencies 
with respect to the DC case and its value is enhanced by 
a factor w A~5 ( see Fig.|5Jl, roughly in agreement with 
experimental datai 

The finite temperature behavior of the dynamical con- 
ductivity is as follows: for T <C CI we get only a small 
correction to the damping term in (Q of order (T/f2) 5 / 3 
while for < T essentially the £7 = result applies 
and the conductivity will be frequency independent and 
equal to the DC value. Therefore upon temperature in- 
crease, the MIC peak is expected to shrink asymmetri- 
cally and eventually disappears from the spectrum, a be- 
haviour consistent with experiments^. The limits of va- 
lidity of the approximations involved in the calculation of 
the spinon correlation functions are m 2 s > c— > m s Q^ 
where the lower bound comes from the effectiveness of 
the saddle point at \x\ ~ Q^ 1 - When expressed in terms 
of temperature, this yields a range between a few tens 
and few hundreds of Kelvin. We expect that the up- 
per limit corresponds to a crossover to a new "strange 
metal" phase, analyzed in 1 ^, where the 7r-flux lattice is 
melt and the "metallic" , linear in T resistivity is recov- 
ered. In some sense the pseudogap phase is on the "insu- 
lating" side of the MIC, and the description adopted here 
is a rather good approximation near the MIC, but the 
"high temperature asymptotics" ~ T 1 / 4 is not correct, 
although it reproduces at lower temperature the inflec- 
tion point in resistivity found experimentally. The above 
calculations do not take into account the holon contri- 
bution to the physical conductivity, but that is of order 
A 4 / 3 , hence negligible for small cutoff A. 

Let us finally discuss how the in-plane resistivity 
anisotropy founc&£ in untwinned single crystals of 
La2_ a; Sr a ;Cu04 (x = 0.02 — 0.04) can fit into our scheme. 
The neutron scattering experiments have revealed in- 
commensurate magnetic structure in lightly doped LSCO 
samples (S < 0.05) Unlike the superconducting LSCO 
compounds where the deviation of the clastic magnetic 
scattering peaks from (jr, 7r) is along the a, b directions 
in the tetragonal basisji^ these peaks are rotated by 



45 degrees around (7r,vr), i.e., they are located along 
the b* axis in the orthorhombic basis. Moreover, from 
the half-width of the scattering peaks one can deter- 
mine the magnetic correlation length in different direc- 
tions. As a big surprise, one finds the correlation length 
strongly anisotropic. In particular, for 5 = 0.024, £' a , = 
94.9A, f£. = 39.9A& The authors of Ref. interpreted 
this result as due to stripe formation along a-axis, but no 
quantitative argument was given. This behavior is fully 
consistent with the magnetic susceptibility anisotropy, 
observed in untwinned lightly doped LSCO crystals (up 
tp 3% doping) >i2i We do not have a quantitative micro- 
scopic theory to consider the anisotropy of the AF corre- 
lation length £, but we can see how such anisotropy can 
be included in our scheme and explore its consequences. 
Suppose the hole distribution is anisotropic (which may 
come from the underlying stripe structure), say the aver- 
age distance between holes is bigger along the a axis, so 
does the distribution of vortices on the AF background. 
To use the nonlinear-cr model treatment of spinons, we 
can rescale the spatial coordinates. The result will be al- 
most the same as in the isotropic case^ except for a coef- 
ficient a in the spinon dispersion ^m 2 . + a 2 v 2 k 2 + v 2 k 2 , 

which reflects the ratio of the AF correlation lengths in 
different directions. Since the spinon mass m s is in- 
versely proportional to the correlation length, we can 
effectively interpret this as m s ^ a < m s ^. To calculate 
the anisotropic conductivity we need to modify the entire 
scheme. However, the major effect can be grasped with- 
out detailed calculation. The diamagnctic susceptibility 
X and Landau damping k due to holons will change, but 
only very slightly, since they come from angular integra- 
tion. On the other hand, in the saddle point calculation 
of the path integral the effect is more pronounced, so the 
values of integrals /, g and hence the numerical factor 
c, mentioned above, will also change. In our spin-charge 
gauge approach the combination cxm 2 = | In 6\r is cru- 
cial in determining the MIC scale. If we assume that the 
basic results of our theory developed for the 2D isotropic 
model survive the generalization to anisotropic case as 
outlined above, one would anticipate the parameter r to 
be also anisotropic. In view of the AF correlation length 
anisotropy we expect r a /rb < 1. As a consequence, the 
peak in a a will be shifted to lower temperature with re- 
spect to (Tf,, as follows from Eq.Q; the anisotropy ratio 
o~ a /o~b will show a sharp increase near the MIC and sat- 
urates as T — > 0, in agreement with experiments^ The 
same phenomenon occurs for the AC conductivity at low 
temperatures, where the factor A makes the anisotropy 
ratio even bigger. To estimate this enhancement, we ex- 
tract the ratio r a /n, by fitting the DC data (the extracted 
value r a /rb = 0.725). We can then use © to evaluate 
the corresponding anisotropy ratio for AC conductivity 
without introducing any additional parameters; a com- 
parison with the experimental curve is shown in Fig. [3] 
This anisotropy is less pronounced than the experimen- 
tally observed value for AF correlation length, as quoted 
above for 8 = 0.024, but this is consistent with the above 
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scheme, where we expect that part of the anisotropy ef- 
fect in the combination cxm 2 s has already been cancelled 
by other effects. 

It is true that these peculiar features in the in-plane 
conductivity are "due to modifications of the dynamics 
of the metallic carriers, and not due to the opening of 
a charge gap" *£ since there is no charge gap in doped 
Mott insulators. What we have shown here is that these 
modifications are due to the presence of a gap in the 
spin excitations and its competition with the dissipation 
which is different from the disorder induced localization. 
For the same reason, the statement of only one mobility 
edge in the scaling theory^ does not apply here. As 
we learned from the authors of Ref. Qr^ the anomalous 
behavior of AC conductivity is observed only up to 6% 
doping. It is understandable that the anisotropy due to 
stripes is not present since they are rotated by 45 degrees 
beyond 6% doping, and their orientation is alternating 
between a and b directions in adjacent layers. However, 
the disappearance of the low-frequency peak cannot be 
explained by the stripe interpretation. On the contrary, 
this is very natural in our interpretation since the MIC is 



not observed in samples beyond 6% doping in the absence 
of magnetic field. Now we make an explicit proposal: To 
do the AC experiment in the presence of a magnetic field 
which would suppress superconductivity and reveal MIC. 
If the deviation from the Drude behavior reappears, that 
would be a confirmation of our interpretation. 

To conclude we have shown that the peculiar in-plane 
anisotropy of DC and AC conductivity observed in the 
lightly doped cuprates can be explained in a unified, 
self-consistent manner within the spin-charge gauge ap- 
proach, and the key ingredient is to attribute the MIC to 
the correlation effect. The anisotropy of the AF correla- 
tion length, and consequently the MIC scale provides a 
rather natural explanation of the observed conductivity 
anisotropy, being most pronounced in the limit Q, T — > 
which is very difficult to explain based only on the stripe 
existence. In fact, this is a crucial experiment to distin- 
guish the disorder- and correlation-induced MIC. 
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for sending us the paper (Ref. 0) prior to publication and 
very helpful correspondence. 
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FIG. 1: Theoretically calculated frequency dependence of the 
AC conductivity for different dopings: 8 = 0.03 (full line), 
5 = 0.04 (dashed) and 5 = 0.05 (dotted). 
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FIG. 2: Calculated frequency dependence of the AC conduc- 
tivity for 5 = 0.03. Also shown is the corresponding DC 
conductivity as a function of temperature (in cm -1 ). 
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FIG. 3: Calculated frequency dependence of the AC conduc- 
tivity anisotropy ratio for S = 0.03. The corresponding DC 
ratio as a function of temperature (in cm -1 ) is shown with 
dashed line. Inset shows the corresponding DC and AC data, 
taken from Refs0and0, respectively. 



